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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


Overview 

Optically-trapped,  strongly-interacting  Fermi  gases  are  models  for  exotic  strongly- 
interacting  systems  in  nature.  For  this  reason,  tabletop  experiments  with  strongly- 
interacting  atomic  Fermi  gases  can  provide  measurements  that  are  relevant  to  all 
strongly-interacting  Fermi  systems,  thus  impacting  theories  in  intellectual  disciplines 
outside  atomic  physics,  including  materials  science  and  condensed  matter  physics 
(superconductivity),  nuclear  physics  (nuclear  matter),  high-energy  physics  (effective 
theories  of  the  strong  interactions),  astrophysics  (compact  stellar  objects),  the  physics 
of  quark-gluon  plasmas  (elliptic  flow),  and  most  recently,  string-theory  (minimum  vis¬ 
cosity  hydrodynamics). 

This  program  has  focused  on  strongly-interacting  Fermi  gases  conhned  in  a  standing- 
wave  CO2  laser  trap.  This  trap  produces  a  periodic  quasi-two-dimensional  pancake 
geometry,  with  conhnement  induced  pairing  between  fermions  of  opposite  spin,  many- 
body  pairing  arising  from  the  Pauli  exclusion  principle,  and  quantum  pressure  arising 
from  the  tight  axial  conhnement.  This  geometry  provides  a  rich  new  paradigm  for 
testing  state  of  the  art  many-body  theories  of  strongly  correlated  systems.  In  addi¬ 
tion  to  this  research,  we  continue  to  be  interested  broadly  in  the  thermodynamic  and 
hydrodynamic  transport  properties  of  strongly  interacting  Fermi  gases. 

Findings 

In  the  past  funding  period,  we  made  a  major  breakthrough  in  our  studies  of  quasi- 
two-dimensional  Fermi  gases,  producing  the  hrst  spin-imbalanced  quasi-2D  Fermi  gas 
and  measuring  its  thermodynamic  properties. 

A.  Spin-Imbalanced  Quasi-Two  Dimensional  Fermi  Gases. 

In  the  experiments,  by  using  a  standing  wave,  as  noted  above,  we  produce  a  pe¬ 
riodic  pancake-shaped  trapping  potentials  with  a  5.3  p,  spacing.  With  200-800  atoms 
per  site,  we  use  site-resolved  imaging  to  study  a  quasi-two-dimensional  geometry. 
For  a  true  two-dimensional  (2D)  system,  the  Fermi  energy  Ep  in  the  weakly  conhn- 
ing  transverse  direction  should  be  small  compared  to  the  harmonic  oscillator  energy 
level  spacing  in  the  tightly  conhning  axial  (z)  direction,  i.e.,  Ep  «  hvz  for 
a  true  2D  system.  Here,  the  atoms  occupy  a  quasi-continuum  of  transverse  states. 
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but  they  are  all  in  the  ground  axial  state.  In  contrast,  when  Ep  »  a  quasi¬ 
continuum  of  axial  as  well  as  transverse  states  is  occupied,  and  the  gas  becomes 
three-dimensional  (3D).  Our  recent  experiments  are  performed  in  the  intermediate 
regime,  where  Ep  ~  hz/^.  In  this  case,  the  system  is  far  from  3D,  but  is  not  quite 
2D.  Enhanced  superfluid/superconducting  transition  temperatures  are  predicted  for 
condensed  matter  systems  in  this  regime,  as  well  as  for  quasi-2D  Fermi  gases,  making 
measurements  in  this  regime  particularly  important  for  comparison  with  predictions. 

Our  previous  work  on  radio-frequency  spectroscopy,  [Zhang  et  al,  Phys.  Rev. 
Lett.  108,  235302  (2012)],  revealed  that  the  spectra  are  not  explained  by  2D-BCS 
theory.  Instead,  a  polaron  model  of  the  two-component  gas  could  explain  the  spectra. 
However,  it  was  unclear  how  the  system  evolved  from  true  2D  behavior,  where  2D- 
BCS  theory  is  valid,  to  behaving  as  a  gas  of  spin-up  and  spin-down  polarons,  where 
each  spin  is  surrounded  by  particle-hole  cloud  of  the  opposite  spin.  We  believed  that 
measuring  the  cloud  prohles  in  a  spin-imbalanced  mixture  would  shed  light  on  this 
problem. 

In  both  the  2D-BCS  and  2D-polaron  regimes,  the  effective  interaction  strength  is 
determined  by  the  2D  conhnement-induced  dimer-pairing  energy  Ef,  and  the  trans¬ 
verse  Fermi  energy  Ep.  Since  Ep  oc  h^/{mL\)  oc  with  L±  the  transverse 

interparticle  spacing  and  Eb  oc  with  Id  the  dimer  size,  the  ratio  Ep/Eb 

sets  the  ratio  of  the  Lj_  to  Id-  For  Ep  «  Eb,  the  dimer  is  small  compared  to  the 
interparticle  spacing.  In  this  case,  we  expect  the  dimers  to  behave  a  small  inde¬ 
pendent  molecules,  without  many-body  effects.  Increasing  Ep/Eb,  the  dimer  size 
becomes  larger  than  the  interparticle  spacing.  In  this  case,  we  expect  to  see  many- 
body  changes  in  the  pairing  energy  and  in  the  thermodynamic  properties.  For  pairing 
in  the  weakly  interacting  regime,  we  expect  2D-BCS  theory  to  be  valid.  However, 
we  hnd  that  2D-BCS  theory  predicts  that  the  pairing  energy  will  be  exactly  Eb  and 
that  the  spatial  prohle  for  a  50-50  mixture  of  two  spin  states  will  be  identical  to 
that  of  an  ideal  gas!  No  many-body  effects  are  predicted  for  measurements  of  the 
radio-frequency  spectra  nor  for  the  spatial  prohles  in  the  2D-BCS  regime. 

In  our  experiments,  we  use  a  mixture  of  the  two  lowest  hyperhne  states  of  fermionic 
®Li.  We  hold  the  number  of  atoms  in  the  majority  state  W  constant,  at  Ni  =  800. 
We  choose  the  number  of  atoms  in  the  minority  state  N2  to  be  0  <  iV2  <  800,  so 
that  0  <  N2/N1  <  1.  Dehning  the  2D-Fermi  energy  of  the  majority,  Ep  =  hv±\/2Ni, 
which  is  held  constant  in  the  experiments,  we  control  the  ratio  Ep/Eb  by  tuning  the 
bias  magnetic  held  near  a  Feshbach  resonance.  We  begin  by  measuring  the  column 
density  of  one  transverse  dimension  x,  as  a  function  of  spin-imbalance  iV2/iVi-  Fig.  1 
shows  the  results  for  two  diherent  bias  magnetic  helds,  at  the  Feshbach  resonance  (832 
G)  and  below  the  Feshbach  resonance  (775  G),  where  the  2D-dimer  pairing  energy  is 
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increased  by  an  order  of  magitude. 
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Figure  1:  Measured  column  density  profiles  in  units  of  Ni/Rtfi  at  832  G,  for  Ep/Eb  = 
6.6  (left  panel)  and  at  775  G,  for  Ep/Eb  =  0.75  (right  panel)  versus  N2/N1.  Green: 
1-Majority;  Red:  2-Minority.  Blue-dashed:  Golumn  density  difference.  Each  prohle  is 
labeled  by  its  N2/N1  range.  For  the  density  difference,  the  flat  center  and  two  peaks 
at  the  edges  are  consistent  with  a  fully  paired  core  of  the  corresponding  2D  density 
prohles.  These  features  are  more  prominent  for  the  higher  interaction  strength  (right 
panel) . 

We  see  in  the  hgure  a  curious  behavior  in  the  column  density  difference,  which 
appears  to  oscillate  for  N2/N1  ~  0.34  for  the  larger  pairing  energy  Eb-  To  exam¬ 
ine  the  behavior  further,  we  measure  the  cut-off  radii  of  the  cloud,  in  units  of  the 
Thomas- Fermi  radius  of  the  majority  Rppi  =  \j‘^Ep/{mu‘l^,  which  is  held  hxed  in 
the  experiments.  When  the  minority  fraction  is  small,  the  cutoff  radius  of  the  major¬ 
ity  density  approaches  Rppi  as  it  should.  However,  as  the  minority  concentration  is 
increased,  both  the  minority  and  majority  spatial  prohles  are  substantially  distorted, 
due  to  the  strong  interactions  between  the  two  components.  Fig.  2. 

We  hnd  that  the  cloud  radii  deviate  strongly  from  predictions  for  an  ideal  Fermi 
gas,  where  Ri/Rppi  =  1  and  R2/RTP1  =  {N2/NiY^^.  Further,  the  spatial  prohles 
predicted  by  2D-BGS  theory  for  a  balanced  mixture  are  identical  to  that  of  an  ideal 
Fermi  gas,  black  circle,  upper  right  in  Fig.  2,  which  is  in  strong  disagreement  with 
the  data. 

As  a  polaron  model  explained  the  radio-frequency  spectra  obtained  in  our  previ¬ 
ous  program  [Zhang  et  al,  Phys.  Rev.  Lett.  108,  235302  (2012)],  we  again  resorted 
to  a  physically-motivated  polaron  model  to  explain  the  cloud  prohles.  At  zero  tem- 
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Figure  2:  Majority  radii  (upper-blue)  and  minority  radii  (lower-red)  in  units  of  the 
Thomas- Fermi  radius  for  the  majority  for  Ep/Eb  =  6.6  (left  panel),  Ep/Eb  =  2.1 
(middle  panel)  and  Ep/Eb  =  0.75  (right  panel).  Dots:  Data;  Dashed  lines:  Ideal 
Fermi  gas  prediction;  Black  circle  upper  right:  2D-BCS  theory  for  a  balanced  mixture; 
Solid  lines:  2D  polaron  model. 


perature,  the  free  energy  density  /  is  equal  to  the  energy  density.  For  an  imbalanced 
mixture,  with  N2  «  Ni,  we  assume  the  2D  energy  density  is 

f  =  ^niepi  +  ^n2ep2  +  n2Ep{2).  (1) 

Here,  the  hrst  two  terms  are  the  energy  density  for  a  noninteracting  gas  and  the  last 
term  is  the  energy  density  for  minority  polarons  in  state  2,  which  arises  from  scattering 
in  the  bath  of  majority  atoms  in  state  1;  ni^2  and  epi^2  are  the  corresponding  2D- 
densities  and  2D-local  Fermi  energies.  The  2-polaron  energy  per  particle  Ep{2)  = 
?/m(?i)eFi,  where  qi  =  epi/Eb-  For  simplicity,  we  use  an  analytic  approximation  for 
ym{qi)  =  — 2/log(l  -|-  2gi),  due  to  Klawunn  and  Ricatti  [Phys.  Rev.  A  84,  033607 
(2011)],  which  interpolates  between  the  Fermi-polaron  and  Bose-dimer  regimes.  From 
Eq.  1,  we  directly  obtain  the  local  chemical  potentials,  =  df/drii  and  /i2  =  Of  /dn2 
and  the  corresponding  local  2D  pressure  p  =  rii  ni  +  n2  —  f  ■  The  local  chemical 
potentials  determine  the  spatial  profiles  in  the  trap,  since  /ii,2(p)  +  f4rap(p)  =  hi, 2(0). 
The  solid  curves  in  Fig.  2  show  the  predictions  of  the  polaron  model,  which  captures 
most  of  the  features  of  the  data. 

We  hnd  that  the  simple  polaron  model  fails  to  predict  a  phase-separation  transi¬ 
tion  that  we  observe  in  the  2D  density  profile  as  N2/N1  is  increased  beyond  a  critical 
value.  For  N2/N1  <  (A^2/A^i)crio  the  cloud  is  a  uniform,  spin-imbalanced  mixture. 
For  N2/N1  >  {N2/Ni)crit,  the  2D-density  profile  changes  to  a  balanced  core  region. 
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where  the  densities  become  equal  near  the  cloud  center,  with  the  excess  majority  spin 
component  expelled  to  the  edges.  A  better  model  of  the  pressure  is  needed  to  explain 
the  transition,  which  provides  an  important  benchmark  for  new  predictions  of  the 
phase  diagram  in  this  interesting  quasi-2D  regime. 

B.  New  Results  on  Strongly  Interacting  Fermi  Gases. 

Recently,  we  have  made  new  measurements  of  the  radio-frequency  spectra  for  a 
balanced  mixture,  as  a  function  of  the  transverse  Fermi  energy.  Here,  we  tune  the  ratio 
Ep/huz  from  <<  1  where  the  system  is  truly  2D  to  ~  1,  where  2D-BCS  theory  breaks 
down.  Our  initial  measurements  show  that  the  2D-BCS  theory  prediction  is  correct  for 
Ep/huz  <<  1  ,  but  that  the  polaron  model  is  needed  to  £t  the  data  for  Ep/lwz  —  1- 
Our  hypothesis  is  that  in  the  strongly  interacting  regime  near  the  Feshbach  resonance, 
atoms  inside  the  dimer  pair  dimension  R,  can  collide  with  each  other,  rendering  a 
Cooper-like  pairing  picture  incorrect.  Instead,  the  atoms  of  each  species  become 
surrounded  by  a  particle-hole  cloud  of  the  other  species,  i.e.,  a  polaron  gas.  These 
new  experiments  are  made  possible  by  an  optical  lattice  system,  constructed  in  part 
using  support  from  DOE,  which  unites  our  ARO  and  DOE  programs. 

We  also  made  a  major  breakthrough  in  our  studies  of  shear  viscosity,  J.  A.  Joseph, 
E.  Elliott,  and  J.  E.  Thomas,  “Shear  viscosity  of  a  unitary  Fermi  gas  near  the  su¬ 
perfluid  phase  transition,”  Phys.  Rev.  Lett.  115,  020401  (2015),  which  was  selected 
as  an  Editor’s  Suggestion.  In  these  experiments,  we  made  the  first  attempts  to  ex¬ 
tract  local  shear  viscosity  from  cloud-averaged  data  and  provide  the  most  precise 
measurements  to  date  of  the  shear  viscosity  from  nearly  the  ground  state  to  the 
high  temperature  regime.  These  experiments  are  spurring  new  theoretical  efforts  to 
understand  hydrodynamic  transport  in  the  strongly  interacting  regime. 

We  made  the  first  study  of  scale  invariance  in  the  hydrodynamic  expansion  of 
a  unitary  Fermi  gas  in  E.  Elliott,  J.  A.  Joseph,  and  J.  E.  Thomas,  “Observation  of 
Conformal  Symmetry  Breaking  and  Scale  Invariance  in  Expanding  Fermi  Gases,” 
Phys.  Rev.  Lett.  112,  040405  (2014).  In  this  work,  we  also  made  the  first  precision 
measurement  of  the  bulk  viscosity,  showing  that  it  is  consistent  with  zero,  as  predicted 
for  a  scale-invariant  system. 
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